Interchain conductivity of coupled Luttinger liquids and organic conductors 
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We reconsider the theory of dc and ac interchain conductivity in quasi-one dimensional systems. 
Our results are in good agreement with the measured c-axis optical conductivity of (TMTSF)2C104 
and suggest that the c-axis dc-conductivity of (TMTSF)2PF6 in the 150K < T < 300K range 
is dominated by precursor effects of Mott localization. The crossover from a Luttinger liquid at 
high energy to a Fermi liquid at low energy is also addressed, within a dynamical mean-field theory. 
Implications for the inter-chain resistivity and Drude weight in the Fermi liquid regime are discussed. 



Inter-chain (or inter-plane) coherence is a key issue for 
the physics of anisotropic strongly correlated systems. It 
is of special importance for the quasi-one dimensional or- 
ganic conductors TMTSF2X and TMTTF2X which are 
systems of weakly coupled chains . Since the one di- 
mensional (ID) interacting electron gas has non Fermi 
liquid (Luttinger liquid) properties one expects in 
these compounds a crossover from Luttinger liquid be- 
havior to a coherent behavior as the temperature or fre- 
quency is lowered. The latter can be either a metallic 
three-dimensional conductor (presumably Fermi-liquid 
like), or a phase with long-range order. 

Although there is general agreement on this qualita- 
tive scenario, there is considerable debate on how pre- 
cisely this crossover takes place As was pointed 
out early on |^ , the interactions renormalize the effective 
interchain hopping of weakly-coupled Luttinger chains. 
This leads to a crossover scale E* oc t(t±/t)^/^^^°'^ = 
tj_{tj_/t)°'^^^~'^^ for the onset of coherence along one of 
the transverse directions with interchain hopping t±. In 
this expression, t is the intra-chain hopping and a = 
\{Kp -\- l/Kp) — I is the Fermi surface exponent. The 
scale E* can be much smaller than the bare interchain 
hopping t±. With this idea in mind, analysis of NMR 
data [p|JlO| led to believe that the ID Luttinger liquid 
(LL) to Fermi liquid (FL) crossover took place in the 
(TMTSF)2X family at a low temperature, of the order of 
10 — 15K. This scale seemed consistent with the above 
estimate of E* since the interchain hopping along the b- 
axis is t}) ^ 300if , an order of magnitude smaller than 
the hopping ta along the chains. Recent analysis of the 
transport and optical properties ||ll[] have led to a re- 
examination of this commonly accepted view. Measure- 
ment of the optical conductivity along the chains agrees 
with a LL picture at high enough frequency and provides 
a determination of the LL exponent Kp |12 1^ (of order 
Kp ~ 0.23). This value is consistent with the temper- 
ature dependence of Pa{T) in the lOOK — 300K range 
(suitably corrected for thermal expansion [Q). However, 
these studies and most notably the recent measurements 
of the interchain dc conductivity jlj] revealed that the 
ID regime only extends down to ^ lOOK, a much higher 



temperature scale than the early estimates based on the 
NMR. Furthermore, it is clear from the optical data that 
another important energy scale is present in these com- 
pounds, associated with a Mott gap signaled by the onset 
of optical absorption at ~ lOOcm"^ ~ 150K. Transverse 
hopping is actually responsible for the metallic character 
of the (TMTSF)2X compounds, which otherwise would 
be Mott insulators [ pl]JT5| . These findings raise a num- 
ber of important and largely unanswered questions such 
as the precise nature of the dimensional crossover, the 
origin of the anomalies observed for 10 K < T < 100 K 
(calling in particular for a reinterpretation of the NMR 
in this range), and the actual range of the LL behavior. 

In this letter, we focus on inter-chain transport and op- 
tical response, which are ideal probes of these issues. We 
first reconsider the theoretical analysis of the w and T- 
dependence of the inter-chain conductivity for the sim- 
plest case of weakly coupled metallic Luttinger liquids. 
We find a power law of the temperature or the frequency 
a±(T,Lu) oc (T, cj)^"^^. The exponent is at variance with 
previous theoretical predictions [g Ui^Ji^ , but is in rea- 
sonable agreement with recent measurements of c-axis 
optical conductivity of the (TMTSF)2C104 compound 
[pTf . We then show that understanding the observed 
T-dependence of c-axis transport in (TMTSF)2PF6 re- 
quires to take into account the physics associated with 
the Mott gap apparent in optical measurements. Finally, 
we consider a dynamical mean-field theory (DMFT) for 
the description of the dimensional crossover, and draw 
some consequences for the inter-chain conductivity in the 
low-temperature FL regime. 

We consider the hamiltonian of coupled chains: 
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in which stands for each individual chain to, and 

the inter-chain hopping t±_ will be first considered, for 
simplicity, to be identical for all transverse directions. 
Using the DMFT approach detailed below, a general for- 
mula (Eq. ^ for the transverse conductivity can be estab- 
lished. When the inter-chain hopping tj_ can be treated 
in perturbation theory to lowest order, this formula re- 
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duces to that estabhshed in Ref. ||TJ] within the tunnehng 
approximation pq]: 



Reai(w,T)f'=''* cx 4 / dx / duj' Aid{x,i^ + 1^') x 



(2) 



where A\y)(x^ll}) is the one-electron spectral function of 
a single chain, and / the Fermi function. Let us consider 
first the simplest case in which each decoupled chain is a 
Luttinger liquid characterized by a given value of Kp , and 
all intra-chain umklapp processes can be neglected. For 
the values of Kp of interest here, the transverse hopping is 
then a relevant perturbation in the renormalization group 
sense, so that it is legitimate to treat it perturbatively 
only at sufficiently high energy or temperature. Eq. (|^) 
thus applies only when kT ^ E* or hu) ^ i?*, where 
E* = t(t_L/t)^/(^~") is the single-particle crossover scale 
mentioned above In this high-temperature/high- 

energy regime, the system can be considered to be quasi 
one-dimensional. Using known properties of the spectral 
functions in a Luttinger liquid, we then obtain from (H): 



Re(Ti(tj,T)P'=''' cx 



t 
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In this expression, is a scaling function such that 



Fq(0) = const and > 1) 



Thus, we find 



Recrj.(T > w) cx T^"-! for kT E* , while at high- 
frequency {fiuj :$> E*), we obtain cr±{LJ T) cx w^"^^. 

This behavior of a±{T) differs from the conclusions 
of previous authors. In the authors used a somewhat 
different approach than the direct Kubo formula and con- 
cluded that aj_{T) cx T^". In ju), it was noted that this 
temperature dependence also results from provided 
that the spatial integral is cutoff at a scale of the order 
of the lattice spacing. According to these authors, this 
cutoff is justified by the incoherent nature of the propa- 
gation of a physical electron within each chain (in which 
the true quasiparticles are spinous and holons). In our 
opinion however, this information is already encoded in 
the spectral function of the physical electron Aid(x,ll!) 
and no extra cutoff has to be introduced in (||). This 
conclusion would of course be changed in the presence 
of specific forms of intra-chain disorder such as a random 
forward scattering, which would spoil momentum conser- 
vation and bring in a natural cutoff of the order of the 
mean-free path. Note however that more realistic forms 
of disorder would be signaled by localisation effects. 

We now turn to a comparison with experiments on 
(TMTSF)2X. It is increasingly recognized that the prop- 
erties of these compounds result from the commensurate 
quarter-filling of the band [Tl|jr^ , |l^ . Transport and opti- 
cal measurements lead to consistent values of the LL pa- 
rameter Kp ~ 0.23. Note that Kp < 0.25 is the condition 
upon which the quarter-filled umklapp process becomes 



relevant. This suggests that if the transverse hopping 
was entirely suppressed in these compounds, they would 
actually be one- dimensional (quarter filled) Mott insu- 
lators. Fortunately, optical experiments extend to high 
frequencies compared to the Mott gap. Thus the com- 
plications associated with the proximity to the Mott in- 
sulating state can be treated perturbatively, making the 
comparison to Eq.(^) relevant. Recently, optical stud- 
ies of (TMTSF)2C104 have been performed 0, and the 
high-frequency optical conductivity CTc(w) along the least 
conducting c-axis was found to be a broad, very slowly 
increasing function of frequency in the range 50cm^^ < 
Lo < lO^cm"^. This slowly increasing continuum is essen- 
tially independent of temperature for lOK < T < 300K, 
and is in good agreement with our result ac{(^) cx uj'^"~^. 
The exponent 2a - 1 = {Kp + l/Kp)/2 - 2 depends 
quite sensitively on the value of Kp {2a — 1 = 0.12 for 
Kp = 0.25 and 2a - 1 = 0.6 for Kp = 0.2). The latter 
value 2a — 1 ~ 0.6 provides a reasonable fit to the ex- 
perimental data of llj]. Note that values of Kp in the 
0.2-0.25 range are quite consistent with the other exper- 
imental findings discussed above. In contrast, the w^" 
law advocated in [|6| Jl^Jl^ would require Kp ~ 0.35 to 
explain these data, a value quite inconsistent with that 
obtained from the longitudinal optical conductivity. Fi- 
nally, an important issue in view of the success of this 
simple theory is whether it can also account for the be- 
havior of the optical conductivity along the b-axis. Since 
th/ta — 0.1, an incoherent regime should also apply there 
at those high frequencies. Preliminary analysis of the 
(7b data of p7[ | seems in agreement with such a behavior, 
but more detailed comparison with the data is needed. 

We now discuss the measurements of the c-axis trans- 
port in (TMTSF)2PF6 ll|,|2l|. Because of the large ther- 
mal expansion, these were performed at several pressures 
so that a volume correction could be made. This yields 
data for the c-axis resistivity Pc{T) which can be inter- 
preted as corresponding to a compound with constant 
lattice parameters. This volume-corrected Pc{T) displays 
a marked increase by more than a factor of 3 upon cool- 
ing from T = 350K down to T = lOOK (followed by a 
rapid downturn at lower temperature which presumably 
marks the onset of inter-chain coherence). Such a large 
variation of Pc{T) cannot be explained by the dependence 
Pc cx iyT2Q!-i fQ^Q(j above with acceptable values of Kp 
(which correspond to rather small values of 2a — 1). In 
||l4| , it was fitted to a 1/T^-^ power law, which happens to 
be in reasonable agreement with the 1/T^" dependence 
advocated by these authors. However, the temperature 
range in which these data are obtained is such that kT is 
never significantly larger than the "would-be" ID Mott 
gap which can be extracted from the optical measure- 
ments (100 — 200cm~^) leading to transport gaps (half 
the optical gap) in the range A ~ 75 — 150K. As a re- 
sult, the proximity to the incipient Mott insulator cannot 
be neglected, and we propose that this effect is the one 
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responsible for the increase of Pc(T). Indeed, we have at- 
tempted a phenomenological fit by an activated behavior 
Pc{T) cx l/T2"-iexpA/T. and found that it provides 
a fit of the data of Ref. [ p"4|j2l]] which is as satisfactory 
as the power-law fit used in |14[|. The value of A used 
in this fit is of the same order of magnitude than the 
Mott gap seen in optical measurements. The discrep- 
ancy between the theoretical result, Eq.(||) and the exper- 
imentally observed behavior of Pc{T) is, in our opinion, 
a strong indication that an explanation based on inco- 
herent tunneling between LL chains in a purely metallic 
regime is insufficient and that the physics of Mott locali- 
sation plays an important role in the temperature range 
150K < T < 300K. 

A more refined theory of the interplay between inter- 
chain hopping and Mott localisation is clearly needed, 
including both the transverse hopping and the quarter- 
filling umklapp scattering |22j. Due to the umklapp, the 
effective Kp will be renormalized towards smaller values 
as temperature is decreased, hence leading to larger val- 
ues of pc, as expected. The main issue is whether the 
intra-chain resistivity is also strongly affected by these 
renormalizations. Indeed, no sign of an upturn in Pa{T) is 
observed experimentally for the TMTSF2X compounds. 
Since pa{T) oc glj^ conceivable that the 

decrease in Kp is almost compensated by an associated 
decrease of the effective 51/4 due to t±. This however, 
remains to be verified and is probably the most chal- 
lenging issue for the picture that we propose. Transport 
in the TMTTF2X family is interesting in this respect. 
These compounds are insulators with a rather large Mott 
gap, and indeed both Pa{T) and Pc{T) display an acti- 
vated behavior at low temperature. Above some tem- 
perature scale however, Pa(T) becomes metallic, while 
Pc{T) still displays insulating behaviour (albeit slower 
than activated). Studies of the transport properties of 
these compounds over a more extended range of pressure 
and temperature (and in particular a detailed compari- 
son of the energy scales appearing in 170(0;), Pa{T) and 
Pc(T)) would certainly prove very interesting. 

Let us finally address the issue of the crossover between 
a Mott-Luttinger liquid at high energy/temperature and 
a coherent (presumably Fermi-liquid) state at low energy. 
The main theoretical difhculty is that this crossover is as- 
sociated with the breakdown of perturbative expansions 
in tj^. On the other hand, perturbative approaches in the 
interaction strength could be used in the FL regime, but 
fail in the LL regime. The necessity of treating t± in a 
non-perturbative manner was recently noted by Arrigoni 
psf , who used a limit of infinite transverse dimensionality 
as a guidance for resumming the perturbation series in 
t±. Here, we show that this limit provides a generalized 
dynamical mean-field theory (DMFT) framework for the 
description of the crossover, and we draw some conse- 
quences for inter-chain conductivity. The limit of infinite 
dimensions has been used extensively in recent years to 



formulate a controlled DMFT of strongly-correlated sys- 
tems 1 24 2^ . We consider each one-dimensional chain to 
be coupled to zj^ — > cx) nearest neighbor ones, with the 
transverse hopping scaled as t± = t±/y/z±, so that the 
non-interacting density of states in the transverse direc- 
tion D{€±) = J2k± '^[^-L ^ ^(''^-l)] remains well-defined. It 
is easily shown that in this limit the self-energy becomes 
independent of transverse momentum fcx: E = !](/(;, iWn) 
{k is the momentum along the chains). Furthermore, 
the calculation of E reduces to the solution of an effec- 
tive one-dimensional problem specified by the effective 
action: 



drdr' c+{T)g^\i -3.T- t')c,,(t') 



- r drHr^[{c,,,ct}] 
Jo 



(4) 



where is the interacting part of the on-chain hamil- 
tonian in (^, while Go is an effective bare propaga- 
tor which is determined from a self-consistency con- 
dition. The latter imposes that the Green's function 
G{i — j,T — t') = —{c{i,T)c'^{j,T'))oS calculated from 
^off coincides with the on-chain Green's function of the 
original problem, with a self-energy E = G^^^ — G^^ . This 
condition reads: 



G{k,iLUn) = / de± 



' iujn + p - (-k - 5](iw„, k) - e± 



(5) 



The DMFT equations (|j5|) fully determine the self- 
energy and Green's function of the coupled chains. Fur- 
thermore, intcr-chain transport properties simplify in 
this approach. Because the current vertex is an odd func- 
tion of the transverse momentum, vertex corrections drop 
out of the inter-chain conductivity, which reads: 



Rea±{uj,T) cx t'i / dej_D{ej_ 



dk 
2^ 



duj'A{€±,k,uj') 



.Aie,,k..^u')f^-'^~^^^'^-^ (6) 

UJ 

Here, A{e^^k^uj) — ilmG'(e^, fc, w) is the single- 
particle spectral function of the coupled chains system. 
Eq.(^ has a much wider range of applicability than its 
small-tj^ limit, Eq.(||), since it holds all the way from the 
LL to the FL regime. 

Solving quantitatively these DMFT equations is still a 
challenging problem however. In the single-site DMFT 
(corresponding to the usual d = 00 limit |^^) the map- 
ping is onto a single- impurity Anderson model with a self- 
consistent bath, and several techniques have been devel- 
oped to handle this problem. Here however, the mapping 
is onto a self-consistent ID- chain, and the problem is an 
order of magnitude more difficult. Nevertheless, several 
conclusions can be drawn from the above equations even 
in the absence of a full solution. We shall consider again 
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the simplest case of metallic ID LL chains, neglecting 
umklapp scattering and Mott physics. In the ID regime, 
the self-energy behaves as S cx From 
(^, it is clear that the inter-chain hopping becomes rel- 
evant when t± > T,, and we recover the crossover scale 
E* cx At energies smaller than E* , the 

DMFT approach leads to Fermi-liquid behavior. This is 
clear from viewing the effective ID model as a set of An- 
derson impurity models coupled together by the non-local 
part of the "bare" propagator Go- Each of the decou- 
pled Anderson model is a Fermi-liquid at low-energy and 
the non-local hybridization is a smooth perturbation on 
the FL ground state. In other words, the physics below 
E* is smoothly connected from small t±/t to large t±/t. 
When tj^ <C t and for energy scales much smaller than 
the bandwidth (of order t), all single-particle quantities 
should obey a scaling behavior in the variables lu — uj/E* , 
X = xE*/t, f = T/E*. In particular tT.{x,LJ,T) = 
EHj:{x,Cj,f), tG{x,u;,T) = {E* /t±)G{x,u,f) with 
S and G universal scaling functions associated with the 
crossover. A low-frequency expansion of the scaling form 
of E in the FL regime yields a finite quasiparticule residue 
Z (X {t±/t)^ cx E*/t± (shown in |Q to have only 
weak k-dependence along the Fermi surface). Z can be 
much smaller than unity, while the effective mass (spe- 
cific heat) ratio m* /m does not become large because 
dTj/du} and d'E/dk scale in the same manner. Using the 
FL form of E in m), we obtain the T = inter-chain 

2 

conductivity as: a±{Lu,T = 0) = ^^^Y^d{(jj) + a^_^^ {u>) with 
WD±V8 oc ZH\/t = t(ti/t)2/(i-"). The total spectral 
weight j dbja = Wp^/S is, from the f-sum rule, pro- 
portional to the inter-chain kinetic energy, and hence to 
i^/t, so that the relative weight carried by the Drude 
peak is small, of order uj'j^^/uj'l,^ cx (t^/O^"^^"^""^ = 
This is in reasonable agreement with the optical data on 
the b-axis reflectivity of (TMTSF)2PF6 for which a 
Drude peak is seen at low temperature, with t^l^^/tjpj^ ~ 
0.03 (while tb/ta — 0.1). Scaling considerations also lead 
to an inter-chain resistivity p±{T)/po = {t/ E*)R{T / E*), 
with R{x > 1) ^ x^^^" , R{x < 1) and po = 

hVyn/ie^a^) This implies that p^/po = A{T/tf in 
the FL regime, with an enhanced value of the coefficient 
A cx (t/t^)T^. Note that p±^ is typically much bigger 
than the Mott limit poj ^ observed experimentally. 

In conclusion, we have reconsidered in this paper the 
theory of interchain conductivity in quasi-lD systems. 
When applied to organic compounds, good agreement is 
found with the frequency dependence of f7c(w). while c- 
axis dc transport appears to be dominated by the proxim- 
ity of the Mott gap A. This calls for additional theoreti- 
cal work in the regime where E* cx i(t^/i)^/(^~") < A < 
t^. It would also be valuable to see whether the regime 
A < i?*, which is simpler to analyze theoretically, can 
be realized experimentally in some family of compounds 
(e.g by applying uniaxial stress to TMTTF2X). 
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